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^ , A system comprising two superconducting thin films connected by a point contact is consid- 

O ' ered. Tlie contact resistance is calculated as a function of temperature and film geometry, and is 

' found to vanish rapidly with temperature, according to a universal, nearly activated form, becom- 

\ ing strictly zero only at zero temperature. At the lowest temperatures, the activation barrier is set 

00 , primarily by the superfluid stiffness in the films, and displays only a weak {i.e. logarithmic) tem- 

■ perature dependence. The Josephson effect is thus destroyed, albeit only weakly, as a consequence 
(— I I of the power-law-correlated superconducting fluctuations present in the films below the Berezinskii- 
Q , Kosterlitz-Thouless transition temperature. The behavior of the resistance is discussed, both in 
^ ' various limiting regimes and as it crosses over between these regimes. Details are presented of a 

minimal model of the films and the contact, and of the calculation of the resistance. A formulation 
in terms of quantum phase-slip events is employed, which is natural and effective in the limit of a 
'i good contact. However, it is also shown to be effective even when the contact is poor and is, indeed, 

' indispensable, as the system always behaves as if it were in the good-contact limit at low enough 

, temperature. A simple mechanical analogy is introduced to provide some heuristic understanding of 

. the nearly-activated temperature dependence of the resistance. Prospects for experimental tests of 

^ ' the predicted behavior are discussed, and numerical estimates relevant to anticipated experimental 

I ' settings are provided. 

^ : 
c ■ 

O ■ I. INTRODUCTION 

o : 
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The Josephson effect^ is one of the most direct probes of superconducting long-range order (LRO). Phase coherence 
- and hence a dissipationless supercurrent - can be maintained across a variety of weak links between two bulk 
three-dimensional (3d) superconductors. While one can imagine the decay of a supercurrent across the link via some 
dissipative process such as a phase slip, this would require the phase to "unwind" deep into the bulk superconductors, 
^vq ' which is prevented by the rigidity of the order parameter. 

In low-dimensional systems, fluctuations are more important and lead to a host of remarkable phenomena. Super- 
conductors are no exception to the rule, and a variety of low-dimensional or quasi-low-dimensional systems where 
, superconducting fluctuations play an important role are the subject of much current interest. For example, several 
experiments have investigated a so-called superconductor-insulator transition in superconducting nanowires.^"^ These 
experiments measured the resistance of thin wires connecting two superconducting electrodes. Theoretical work on 
the subject^"^ has explored the possibility of universal length- and width-dependent behavior of such devices. Also, 
C [ superconducting phase fluctuations play an important role in some strongly correlated and disordered systems. Note- 
. worthy examples include the pseudogap regime of the cuprate superconductors^ and disordered two-dimensional (2d) 
' films near a superconductor- insulator transition. 
O ■ A striking aspect of low-dimensional systems is that the fluctuations can destroy true long-range phase order, but 
^ ■ sometimes only weakly so. Rather than a resistive state with only short-range superconducting correlations, one finds 
^ , power-law superconducting fiuctuations, i.e. quasi-long-range order (QLRO), and vanishing bulk resistivity. This is 
JJ"! ■ the case in one-dimensional (Id) superconducting wires at zero temperature, and in 2d films below the Berezinskii- 
rN ' Kosterlitz-Thouless (BKT) phase transition at temperature Tbkt- 

What is the fate of the Josephson effect for such low-dimensional superconductors? The one-dimensional case 
■ ■ '■ of this question, together with a host of related issues, has been the source of much experimental and theoretical 
interest. ^^^^^ While this work has been pursued in the context of various physical systems, the common thread is that 
the Id conducting ground state is a Luttinger liquid with power-law superconducting correlations and zero resistivity. 
The simplest system consists of two Id wires joined by a point contact. There, provided the superfluid stiffness of the 
wires is greater than a certain critical value, the tunneling resistance vanishes as a power-law in temperature, with 
a universal exponent depending only on the stiffness. This power law is a reflection of the critical nature of the 
zero-temperature Luttinger-liquid ground state of the wires. 

Perhaps the simplest system for studying Josephson physics in two-dimensional superconductors consists of two 
superconducting thin films connected by a point contact, as illustrated schematically in Fig. lA. The point contact 



FIG. 1: We consider both wedge (A) and bilayer (B) geometries. The shaded regions are the superconducting films, and the 
dark line joining them represents the point contact. Note the definition of the opening angle 6 in the wedge geometry. 



can be any kind of weak link much smaller than the films themselves. Systems of this kind, where two MoGe films are 
joined by a thin nanowire^'^ or a single narrow constriction,^^ have been the subject of recent transport experiments. 
A related situation has been studied in the experiment of Ref. 23, where Pb films were probed by the tunneling of 
Cooper pairs from a superconducting STM tip. Aside from the work of Kim and Wen,^^ which we discuss further 
below, very little theoretical attention has been paid to this problem. 

Motivated by these experiments, and in anticipation of further work on similar systems, we have studied theoretically 
the d.c. transport across a point contact between two thin-film superconductors. A nontechnical discussion of our 
main results can be found in Ref. 25. We focus on linear-response properties; for example, we consider only the 
resistance given by the voltage response to a very weak current bias. The point contact is characterized by the 
Josephson coupling energy J, or, equivalently, the critical current Ic — 2eJ/h. The films have superfluid stiffness 
Ks = h^ns/m, where Ug is the number of Cooper pairs per unit area of the film, and m is the pair mass. The problem 
can be approached from the opposite limits of a poor contact {J /Kg small) or a good contact [J /Kg large). The 
transport in the poor-contact limit is naturally viewed in terms of Cooper pair hopping events, where a single pair 
hops across the contact from one film to the other with amplitude J. This limit was analyzed by Kim and Wen,^** 
who calculated the d.c. conductance Q{T) as a function of temperature, and found that it apparently diverged below 
a characteristic temperature T* .^^ 

As we argue below, the result of Ref. 24 is in error (except at zero temperature) , and is symptomatic of a breakdown 
of perturbation theory rather than a physical zero-resistance state. We arrive at this conclusion by working in the 
good-contact limit; in fact, one of our main results is that, at low enough temperature, the leading-order perturbative 
calculation in this limit is always correct, even when J <C Kg- Transport in the good-contact limit can be viewed 
in terms of quantum phase slips; these are events where the order- parameter phase winds by 2tt across the contact. 
Phase slips correspond to vortex-hopping processes, in which a vortex enters one of the films near the contact and then 
moves to the other side of the contact, where it leaves the film (see Fig. 2). A different, but qualitatively identical, 
phase-slip process involves the phase at one end of the contact changing continuously in time from to ±27r; this can 
be visualized as a vortex hopping directly across the contact (Fig. 2). The disturbance created by one of these events 
can heal locally, via the reverse process where a vortex hops in the opposite direction. Alternatively, it can propagate 
out into the films, causing a voltage spike across the system of films plus contact, according to the Josephson relation 

where tp is the phase of the superconducting order parameter, and AV and Aip schematically represent the difference 
in voltage and ip, respectively, between two points. Formally, it is possible to work directly in terms of phase-slip 
events, via a duality transformation. One can then calculate the resistance perturbatively in ty, the amplitude for 
phase-slip events to occur. 

We note that the phase slips are quantum in nature because t,y approaches a constant as the temperature is 
lowered to zero; in this limit phase slips are clearly quantum tunneling events. Furthermore, we expect only a 
weak temperature variation in ty as long as T <C Tbkt- However, the physics determining the low-temperature 
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FIG. 2; Depiction of phase-slip processes. In one, a vortex (ellipse with arrow) enters the left-hand superconducting film below 
the contact (dark line with "X" through it), follows the dashed- line path, and exits the film on the other side of the contact. In 
the second, the vortex moves directly through the contact, along the dash-dot line. This process corresponds to a continuous 
change of the phase at one end of the contact from zero to ±27r, while the phase at the other end remains fixed. 



behavior of the resistance is that of classical phase fluctuations within the superconducting films, except at strictly 
zero temperature. The disturbance caused by a vortex hopping across the contact is an excitation of the gapless 
plasmon mode in the films, and it is the physics of the plasmons that determines whether a voltage spike will register 
across the; ciiitirc system. Our detailed considerations show that the plasmon excitations important in determining 
the low-temperature resistance have frequency much lower than kBT/H, and are thus thermally populated (and hence 
in a classical regime) at any nonzero temperature. 

We find that the resistance vanishes at T = 0, the only temperature where there is superconducting LRO in the 
films. For all T > 0, on the other hand, the resistance is nonzero but vanishes following a universal form as T 0. 
This implies that the perturbative poor-contact result is incorrect; at any fixed temperature, the resistance cannot 
decrease as the contact is weakened, so, as the resistance is nonzero in the good contact limit, a vanishing resistance 
is out of the question. The low-temperature form is very nearly an activated exponential, 

where Rq = h/Ae^ is the quantum of resistance, and Ea{T) is a temperature-dependent effective activation barrier: 

Ea (T) = cKs in(;,^j2fcBT) + (2cX,/7r2 J) " 

Here, J « J is an effective Josephson coupling that need not be distinguished from J for most purposes, and Uc 
is a high-energy cutoff set by an appropriate short-distance cutoff in the films; these parameters are discussed in 
more detail below. The superfluid stiffness plays a crucial role in setting Ea{T). The constant c is a universal 
dimensionless quantity given by 

where 9 is the "opening angle" at the point contact, as shown in Fig. lA. For Coulomb interactions between electrons 

in the films we have a = 2, whereas a = 1 if the interactions are screened down to a short-range form. The latter case 
can be achieved in the presence of a superconducting ground plane separated from the films by an insulating layer 
[see Sec. IIIB]. 

As written, Eq. (2) is expected to be asymptotically exact as T ^ 0, in the sense that [In i?(T)]/ [In i?nicasurcd] ^ 1 
as T ^ 0. Furthermore, this formula should be a good approximation whenever EA(T)/kBT > 1. It is important 
to note that the ratio of the resistances themselves (as opposed to the ratio of their logarithms) will not approach 
unity as T ^ 0. This is due to additive corrections to Ea{T) not captured in Eq. (3); the leading such correction 
is proportional to 1 /T In^ T. These corrections, which are expected to depend on J, are recovered if the integrals 
leading to Eq. (2) arc evaluated exactly. This point is discussed further in Sees. VB and VIII, where the integrals are 
written down in closed form. 

Focusing on the barrier Ea{T), the dependence on J drops out at low temperature, where we find the universal 
form: 

cK 
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This form is universal in the sense that it does not depend on the details of the point contact; instead, the important 
temperature scale is set by Kg, which is a property of the films themselves and can be determined independently by 
an inductance measurement.^^ The same result holds in the limit of a good contact (J ^ Kg) for all T <C Ks/ks- 
At the lowest temperatures, then, the system always behaves as if it were in the good-contact limit. 

On the other hand, when J <^ Kg, provided the temperature is not too low, we find that the barrier has the form 



More precisely, for this to hold we require that \n{fkUc/kBT) <C ln{hwc/kBTj), where Tj is a crossover temperature 
defined by 



It should be noted that the condition to be in this limit involves the logarithm of the temperature, rather than the 
temperature itself. In the opposite limit (i.e. low temperature), the resistance always crosses over to the universal 
form of Eq. (5). As shown in Appendix C 2, the form Eq. (6) is expected to be asymptotically exact when the above 
conditions (on J/ Kg and the logarithm of the temperature) are satisfied, provided we also require that J 3> fc^T. 

Therefore, Eq. (3) is an approximate interpolation between two limits where it is exact. However, the order-t^ 
contribution to the resistance is expected to provide an essentially exact result provided ksT <C inm{Kg, J). Our 
result for Ea{T) comes from an approximate evaluation of this contribution, which is exact in the limits discussed 
above. We provide a closed-form expression for the exact 0{tl) contribution in Sec. VB, which can be evaluated 
numerically. The reason that terms of higher-order in ty are not expected to be important is that these terms 
involve larger activation barriers, owing to the strong binding of phase slips in an imaginary-time functional-integral 
formulation; therefore, the 0{ty) term will dominate whenever ksT is smaller than Kg and J, as these energy scales 
set the barrier height. 

The nearly-activated form of Eq. (2) is suggestive of the experimental results of Ref. 22. Furthermore, as discussed 
in Sec. VII, it should be possible to carry out future experiments in a regime where a detailed comparison with our 
theory is appropriate. 



We continue in Sec. Ill with a description of our model, which is a quantum phase Hamiltonian for the films that 
focuses on quantum and thermal fluctuations of the Cooper pair phase field. The contact is modeled by a point-like 
Josephson coupling having strength J. The degrees of freedom in the films are integrated out to obtain an effective 
action involving only the phase difference across the contact. We also give some general arguments in favor of a 
treatment starting from the good-contact limit. 

In Sec. Ill A we discuss some subtleties in treating Coulomb interactions in the wedge geometry of Fig lA, and in 
Sec. IIIB wc discuss a setup where the films and contact are placed in proximity to a superconducting ground plane, 
in order to screen the Coulomb interaction down to a short-range form. In Sec. IV we review the calculation of the 
conductance in the poor contact limit {i.e. working perturbatively in powers of J). 

In Sec. V A we describe the duality transformation used to access the good-contact limit. In Sec. VB, the resistance 
is calculated in this limit, to leading order in ty {i.e. the amplitude for quantum phase slips). Some further intuition 
for the physics behind the resistance formula is provided by a mechanical analogy, which we discuss in Sec. VI. 

In Sec. VII we discuss the prospects for experimental tests of our theory via measurements of point-contact tunneling 
between two superconducting films. We give necessary conditions for such systems to be in a regime where our theory 
applies. We conclude in Sec. VIII with a discussion of some of the open questions raised by our results for a variety 
of different systems. 

Several appendices contain the more technical aspects of the work. Appendix A contains a derivation of the 
imaginary-time action for the point contact, which is obtained by integrating out the low-energy degrees of freedom 
in the films. The analytic continuation used to obtain the necessary correlation functions in real time is discussed 
in Appendix B. The details of the evaluation of the integrals used to obtain the resistance formula of Eq. (2) are 
contained in Appendix C. Appendix D presents the details of a variational approach to the point-contact problem, 
which provides another justification for the good-contact treatment of Sec. V. 
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III. MODEL OF SUPERCONDUCTING THIN FILMS CONNECTED BY A POINT CONTACT 



We consider a system of two superconducting thin films connected by a point contact. Essentially, this can be 
any weak link much smaller than the dimensions of the films. In some cases, especially, for example, if the link is a 
superconducting wire, it is necessary to be sure that one is in a regime where any internal degrees of freedom of the 
link play no role and it can be treated as a zero-dimensional object. We use spatial coordinates r defined separately 
for each film, so that the contact is at the origin r = 0. We are interested in the "wedge" geometry of Fig. lA, which 
is specified by the opening angle 0. Although we also consider the bilayer geometry of Fig. IB, our treatment of that 
case is intended only to illustrate features of the wedge setup in a technically simpler context. In particular, we do 
not include the interlayer Coulomb interaction, which would be important in a real bilayer system. However, it should 
be noted that a real bilayer system can be analyzed along the lines of Sec. IIIB, where we consider a system of films 
in the wedge geometry, separated from a siiperconducting ground plane by an insulating barricir.^^ 

Note that in this Section, as well as Sees. IV through VI and the appendices, we shall set = /cb = 1. 

We are interested in T < Tbkt, which we assume to be substantially below the quasiparticle energy gap. We 
can therefore focus on the Cooper-pair degrees of freedom in the superconducting films. Within the superconducting 
state, the physics is encapsulated in a quantum phase Hamiltonian written in terms of the Cooper-pair phase (pi{r) 
and density ni{r), which satisfy the canonical commutation relations [ipi{r),nj{r')] = i5ijS{i — r'). Here, i,j = 1,2 
label the two films. In the wedge geometry we must impose the condition that supercurrent cannot fiow across the 
film edge; mathematically, this gives the boundary condition V(pi • n = 0, where n is normal to the edge. Due to this 
boundary condition, the case 6 = 2it does not correspond exactly to the bilayer geometry; however, we find identical 
results in these two cases. 

The Hamiltonian for the system is 

n = ni + H2-Jcos{^i{o)~Moj), (8) 

where the last term is the Josephson coupling at the contact, and Hi, 2 are the Hamiltonians for the thin films. These 
take the form 



n, 



= ^J dM'^^ir + Vi, (9) 



where Kg = Ks{T) is the superfluid stiffness in units of energy, and Vi is the interaction potential. For Coulomb 
interactions this takes the usual form 



2 

We also consider the case of short-range interactions, which take the simple form: 

_ 

2K 



[ d^rd^r'^^P^. (10) 
J \r-r'\ 



= ^ / d'r{niir)y, (11) 



where Vg is the superfluid velocity. The short-range case directly describes a two-dimensional superfluid of charge- 
neutral bosons. However, it can also be realized for a superconducting film if the Coulomb interaction is screened 
down to a short-range form. This can be achieved in the presence of a superconducting ground plane, as discussed in 
Sec. IIIB. In both cases we neglect inter-plane density-density interactions, which are expected to be unimportant in 
the wedge geometry. Further justification for this expectation is given below. 

The low-energy excitations in the films are the superconducting plasmons, whose thermal excitation is responsible 
for the QLRO occurring for < T < Tbkt- For screened interactions (and a translation- invariant film), one finds an 
acoustic plasmon having the linear dispersion ui = Vslk]. In the case of Coulomb interactions the plasmon has the 
dispersion lo = 2ey/2nKs\k\. 

Note that we have completely neglected vortex fluctuations within the films. We can justify this by estimating the 
typical number of vortices in a film of linear dimension L. The energy cost of a single free vortex in a 2d superconductor 
is on the order of irKg ln(A/^) plus a core energy -Ecorc- Here, ^ is the coherence length and A = /A-kc^Ks is the 
magnetic screening length. A can be quite large; for example, if Tbkt = 5 K, A(Tbkt) ~ 0.8 cm. Now, the number of 
vortices in the whole film is approximately given by 

JT) (12) 
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At T = Tbkt, the exponent irKg/T takes on the universal value of 2, and, neglecting -Bcoro, Ny « (_L/A)^. Thus, 
for L = 1 cm films and A = 0.8 cm, iV„ is approximately unity at Ibkt and decreases rapidly with temperature. 
Therefore, the typical separation between free vortices in the films is on the order of or greater than the system size 
L, and vortices can indeed be neglected. 

As the modes of the phase field in the films are governed by a harmonic theory, we can integrate them out to obtain 
an imaginary-time action for the phase difference across the contact, (/)(t) = (pi{r = 0,t) — (/32(r = 0,r). Details 
of this procedure, which is straightforward except for the case of Coulomb interactions in the wedge geometry (see 
Sec. Ill A), are given in Appendix A. The result is Z = J T>(j)ex.p{—So — Sj), where 



_ 2cKsT ^ 



\4>{uJn)f. (13) 



Here c is defined in Eq. (4), a;„ e 27rTZ, and we have defined the Fourier transform via 



<l>{uJn) = / rfTe^'^"-0(T) (14) 







cj>{r) = T^e-'^-^4>{oJn), (15) 

where [3 = 1/T. Finally, lOc is a high-frequency cutoff set by the short-distance cutoff in the films, which we take to 
be on the order of the coherence length ^. For example, in the bilayer geometry Wc = t^(|fe| = ^w/ac)- The remaining 
term in the action is simply the Josephson coupling across the contact: 



-J 



/ dT cos(</)(t)). (16) 



We note that the summand of depends only weakly on a;„; apart from the important logarithmic frequency 
dependence, it is simply a |0(w„)p "mass term" for the field. Such a term is obtained for a point contact between 
3d superconducting electrodes, and is thus expected to lead to zero resistance even at finite temperature. It is 
therefore apparent from the form of that the system of films plus contact is in a sense very close to exhibiting 
true superconducting behavior (i.e. a d.c. Josephson effect at finite temperature). It is not, however, clear a priori 
whether the Josephson effect survives intact at finite temperature, or whether it is weakly destroyed. Our analysis 
shows that the latter scenario is the correct one. 

It is instructive to consider the physics at zero temperature, where we argue there will be no quantum phase 
transition as a function of J or other parameters. This can be seen by considering phase-slip instanton events in 0, 
where (j) (f> ±2?: over some localized region in imaginary time. It is useful to recall a different model, where 5*0 is 
replaced by S'q oc J2uj \^n\\<p\'^- In that case, very similarly to vortices in the 2d XY model, phase slips interact via 
a potential proportional to the logarithm of their separation, and there is a zero-temperature Chakravarty-Schmid 
phase transition as the coefficient of S'q is varied. ^^'^^ If the potential is stronger than a logarithm at asymptotically 
large times, the potential energy will always dominate over the entropy, and phase slips will be bound into "neutral" 
pairs {i.e. pairs of one <^ — > -|- 27r phase slip and one 4> ^ (j) — 2tt phase slip). This is exactly what happens in the 
present case, where we find a linearly confining potential (up to important logarithmic corrections). Therefore, the 
system is always superconducting at T = 0. 

The very strong binding of phase slips already indicates that one should be cautious about the expansion in the 
poor contact limit. There will be many events in imaginary time in which a Cooper pair tunnels across the contact, 
but the poor-contact expansion is most reliable when we can think of a dilute gas of such events. On the other hand, 
the good-contact expansion, in which one thinks instead of a dilute gas of phase slips, should be very reliable. 

Some comments are now in order regarding our neglect of density-density interactions between the two films. For 
films of finite size, inter-film interactions will lead to a Coulomb blockade below a temperature that vanishes as the 
film dimensions are taken to infinity. This point is discussed further in Sec. VII; here we concentrate on the limit 
of infinite films. For short-range interactions in the wedge geometry, such interactions are certainly unimportant. 
They will lead to an additive contribution to So proportional to w^l^p. At low frequency this term is dominated 
by the logarithmic term of Eq. (13) that came from the power-law phase fluctuations in the films, and there will 
be no effect on the low-temperature physics. The case of inter-film Coulomb interactions is more subtle. However, 
following Sec. IIIB, we have analyzed the bilayer geometry for films of equal stiffness, including the inter-film Coulomb 
interaction.^^ This is a useful case to consider, because inter-film Coulomb interactions should have a stronger effect 
than in the wedge geometry. It is convenient to exploit the symmetry under interchange of the films and make the 
change of variables (p±{r) = {tfiir) ± ip2{r))/V2, and similarly for the densities. The Hamiltonian decouples into a 
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sum of two terms, each depending only on cither symmetric or antisymmetric fields. The Josephson term couples only 
to the antisymmetric sector, so the symmetric sector docs not contribute to the charge transport. The main effect 
of the inter-film Coulomb interactions is to screen the interaction in the antisymmetric sector down to a short-range 
form, leading to a linearly-dispersing acoustic plasmon. This is quite natural, as the density n_(r) is that of charge- 
neutral dipoles of positively charged Cooper pairs in film one and their negatively charged counterparts in film two. 
The qualitative structure of the effective action is completely unaffected, and, not surprisingly, one now finds a = 1 
as appropriate for a model with only short-range interactions. The effects of inter-film Coulomb interactions in the 
wedge geometry should be even less. Indeed, we expect a = 2 as if the inter-film interactions were absent. This is 
because the special structure of the bilayer, which led via screening to a linearly-dispersing plasmon mode, is absent 
in this case, and the plasmons in both films should still follow u> oc ^/\k\ dispersions. 



A. Coulomb interactions in wedge geometry 



In the bilayer geometry it is completely straightforward to obtain the action Eq. (13) by working in Fourier space. 
The problem is also simple for short-range interactions in the wedge geometry; the kinetic term in Tii is diagonal in a 
basis of eigenfunctions of the Laplacian, and in this case Vi is diagonal in any basis. Coulomb interactions are more 
difficult to treat, as it is not straightforward to simultaneously diagonalize the kinetic and potential terms. We believe 
it is possible to avoid this difficulty without affecting the result by replacing the Coulomb potential by a different 
form V{r,r') that has some of the same properties, but is more amenable to analytic treatment. 

In an infinite 2d plane, the Fourier transform of the Coulomb potential Vc{r) = 4e^/|r| is 

Vc{k) = J (fre-^'-^Vcir) = (17) 

We can write 

Vc = 8^^e\^/^)-\ (18) 
where we define the square root of the Laplacian in a plane wave basis, that is 



V-V2e*'="' = |fe|e*'="'. (19) 

We shall construct a potential in the wedge geometry for which the analog of Eq. (18) is true. 

In the wedge geometry it is natural to work in a basis of eigenfunctions of the Laplacian that respect the V<^ • n = 
boundary condition. In plane polar coordinates {r,ip), this condition takes the form 

d^xir, V = 0) = d^xir, i^ = 6)=Q. (20) 

The eigenfunctions are labeled by a non- negative real number A and a non- negative integer n, and satisfy — V^XAn = 
}?X\n- They take the form: 



j^„/,(Ar)cos(^), n>0, 

X\n{r) = S V , , (21) 
fpoiXr), n = 0. 

Here, Ji,{x) is a Bessel function of the first kind. These eigenfunctions are orthonormal, 

POO p6 

/ drr I dil) x\n{r)x\' 5(A-A'), (22) 

JQ Jo 



and satisfy a closure relation: 



poo 

/ dXj2xxn{r)xxnir') = dir-r'). 



(23) 



With this machinery at our disposal, we define \/— by: 



V^XXn = Ax An- (24) 



FIG. 3: Illustration of the system of superconducting films (white regions) and point contact, separated by an insulating layer 
of thickness d from a superconducting ground plane (shaded region). For the geometry shown, 6 = -k. 



This then leads to the potential: 

y(r,r') =87reM dX^, ''^ , (25) 



which satisfies ^/^-V^V{r, r') = 8ne'^6{r — r'). This function depends separately on r and r', and not only on |r — r'|. 
As the wedge geometry already lacks global translation and rotation symmetries, for imiversal properties it should be 
harmless to sacrifice these in the form of the potential. Furthermore, V{r, r') has the same scaling behavior as the 
Coulomb potential: 

V{sr,sr') = -V{r,r'). (26) 
s 

This can be shown by working with the definition Eq. (25) and observing that 

XA,n(sr-) = -^Xs^nir)- (27) 
Upon changing variables Eq. (26) follows immediately. 



B. Superconducting ground plane 

Here we consider a system of films and point contact in the presence of a superconducting ground plane. As shown 
in Fig. 3, the ground plane is separated from the films by an insulating layer of thickness d. The motivation to 
consider this setup is that the Coulomb interaction is screened with the ground plane in place, which allows a = 1 to 
be realized in a system of superconducting (as opposed to superfluid) films. Furthermore, this setup is quite favorable 
for the design of experiments that can be compared to our theory over a wide range of temperatures (see Sec. VII). 

We model the ground plane as a superconducting film having stiffness Kf ^ Kg, as is appropriate for a thick layer 
of superconducting material. Furthermore, we assume that the insulating barrier is sufficiently high and thick that 
Josephson tunneling between the films and the ground plane can be neglected. For simplicity, we set the dielectric 
constant of the insulating medium to unity (i.e. vacuum). The phase and number fields of the ground plane are 
denoted by ^g{r) and ng{r), respectively. 

The goal is to derive an action like Eq. (13) for the phase difference across the contact. It will turn out to be good 
enough to consider a problem in a simpler geometry, where we have only a single translation-invariant film and the 
ground plane {i.e. no contact). In the limit A'f ^ Kg, we conclude that the geometry of Fig. 3 is well described by a 
model having no ground plane, and the short-range interaction of Eq. (11), with superfluid velocity Vg = y/l&Tre^ Kgd. 

The Hamiltonian for the coupled system of a single film and the ground plane is Wbiiayer = 'Hk + 'Hmt, where 

= f / (V^)^ + ^ / rf^r- (VvP,)^ (28) 

and 

Here we have defined the three-dimensional charge density 



N{R) = 5{z - d/2)n{r) + 6{z + d/2)ng{r), 



(30) 
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where R = (r, z) is the 3d position. In momentum space, 

where it = (fe, kz)- The fcx-integral can be done via contour integration, yielding 

f £k Ane\ ( 1 e-'^l'^l \ / n(-fe) \ 

Wint simplifies in the basis n± = (l/\/2)(n ± n^), where one obtains 

Wint = / ^^[(l + e-''l'=l)|n+(fe)p + (l-e-''l'=l)|n_(fe)p]. (33) 

At this point, it is straightforward to reformulate the problem in terms of an imaginary-time functional integral, 
integrating out n± (fe) to obtain the following effective action for the phase fields: 

Sesm = ^ E / '^r.)M,s{k, iVnM-k, -u;„), (34) 

where $'^(fc,w„) = (<^(fc, a;„), a;„)) and 

\ 167re2sinh(d|fe|) 87re2(l-e-2<il'=l) + -"'s'^ / 

The dispersion of the plasmon modes is given by solving det Mes{k,itj) = 0. This results in a quadratic equation 
for w^, which yields the following two solutions when d\k\ ^ 1: 

(j{k) = Vs\k\, (36) 
ujg{k) = ^J32^^e^{K, + K|)\k\. (37) 

The first of these is a linearly dispersing acoustic mode with velocity 



' \/\<a-Ke^Ksd. (38) 



An analysis of the zero eigenvectors of Meff(fe, iuj{k)) shows that this mode propagates entirely in the film in the limit 
Ks ^ K^. Similarly, the other plasmon mode propagates entirely in the ground plane in this limit. In this limit 
Vs = y/ldne'^Ksd. 

Therefore, when Kg <C K§, the phase correlations in the film will be controlled entirely by the acoustic mode, and 
a treatment starting from a model with short-range interactions should give the correct action for the phase difference 

across the contact in the geometry of Fig. 3. This expectation is substantiated by an explicit calculation, in the 
simpler translation-invariant geometry discussed above, of the correlation function 

C{oJn) = / dT e"^"^ {ip{r = 0, r)(p(r = 0, 0)) . (39) 
Jo 

In general, this correlator for the phase at the point contact completely determines the desired effective action. For 
small Matsubara frequency, the result is 



which reduces to the result for a single film with short-range interactions when K§ Kg 
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IV. POOR CONTACT 



We review the treatment in the poor-contact Umit {i.e. J <C Kg), keeping in mind that such an approach should 

be viewed with caution. Wc can use a Kubo formula to calculate G{T), the linear-response d.c. conductance across 
the contact, to lowest order in J. We shall show, below, that the resulting conductance appears to diverge below 
a characteristic temperature T*. These results agree with those of Ref. 24. This calculation alone is not enough to 
say whether the divergence represents a true zero-resistance state, or simply a breakdown of perturbation theory in 
J. The considerations in Sec. V below demonstrate that the conductance is finite for all T > 0, and, therefore, that 
weak-coupling perturbation theory fails below T*. 

The operator for the current flowing across the contact is 

i = 2eJsm{4>). (41) 

Standard techniques of linear-response theory can be used to derive the Kubo formula relating the conductance to 
the retarded Green's function of the current: 



g{T,u;) = -Re -G^(a;) , (42) 

-OJ 

where 

G'j,{io) = -i[ dte{t)e-'^*{[i{t),i{0)])j=o- (43) 

J —CO 



We shall be primarily interested in the d.c. {u — > 0) limit. 
Straightforward manipulations lead to the expression 



Gjj(u;) = -z(2eV2) / dte 



exp 



(-(<^2-0(O<^(O)))-(i^-i) 



(44) 



Here, we have dropped terms involving exp(— (^^ -|- ^(f)^(O))), which vanish in the thermodynamic limit. Inserting 
the result for the phase correlation function, obtained by analytic continuation in Appendix B 2, we find 

poo 

G^(a;) = -4e2 J2 / dte''''* sin (S(t)) exp (Ait)) , (45) 
Jo 

where 

A{t) = -^J^ °da;(cos(a:t)-l)coth(^), (46) 

In the second equality of Eq. (47) we have dropped a term oscillating at the cutofi^ frequency Wc, which is an unphysical 
artifact of the hard cutoff. We now have, in the d.c. limit, 



/■oo 2 , 

g{T) = Ae^J^ dt t sin (^^) exp {A{t)) . 



(48) 



The integral of Eq. (48) can only diverge due to the contribution from long times. To extract this behavior we 
first need to understand the behavior of A{t) as t — » oo. In this limit, we have 1/t <C T <C Wc, and the dominant 
contribution to A{t) comes from the region of the x-integral where 1/t <^ x <^ T. There, the cosine is oscillatory and 
can be dropped, and coth(/3a;/2) « 2/(3x. We then find that 

A{t)^~^HtT). (49) 
At long times, the sine in Eq. (48) can be expanded to leading order, and the long-time contribution is thus 

_2 2 t2 poo 

g{T)^^^-L dt{tT)--^'^'^''% (50) 
civ,, ./*„ 
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where to is arbitrary. It is clear from this expression that G{T) apparently diverges for 

T<T* = ^. (51) 

We shall show, below, that this divergence is due to a breakdown of perturbation theory, and the conductance is 
in fact always finite for T > 0. 



V. GOOD CONTACT 



Given the apparent divergence of the conductance below T* in the limit of a poor contact, it will be fruitful to 
consider the opposite case of a good contact. Even without knowing the poor-contact result, we should expect an 
expansion about the good contact limit to yield correct results, as discussed in the Introduction and Sec. III. Moreover, 
a variational calculation discussed in Appendix D provides further support that the good-contact limit is always the 
correct starting point at low temperature. In this Section, we first implement a duality transformation that allows us 
to work directly with the phase-slip degrees of freedom. Following that, we use the dual formulation to calculate the 
resistance of the point contact. 



A. Duality transformation 



The starting point for the duality transformation is the partition function for the phase difference across the 
contact, Z = J 2?0exp(— 5*0 — Sj), where and Sj are defined in Eqs. (13) and (16). We proceed by replacing the 
cosine of Sj with a Villain function: 



exp( J cos (/)) ^ ^ exp - ^ (0 - 2nr] f^ . 



(52) 



This is merely a different 27r-periodic potential, so the universal physics should be unaffected. To make this replacement 
in the partition fimction, we need only make r] a function of imaginary time (which we can take to satisfy periodic 
boundary conditions): 



Zv = 



P0 E 

77(t)=- 



exp 



2cK,T 



El0K)l^[ln(^) 



(53) 



We have replaced the Josephson coupling J by the energy J. The new partition function, with the cosine replaced 

by the Villain function, can be thought of as an effective description of the original one, and J as an effective J. 
However, it is better to think of J and J as parameters arising in two (slightly) different effective descriptions of the 
same underlying system. Although we expect that J « J, it is not important to have a more precise relationship 
between these parameters, or even to distinguish them, imless one of them can be directly measured. Indeed, in the 
absence of such a measurement, J must be treated as a fitting parameter in the resistance formula of Eq. (2). 

To specify Zy more precisely, it is useful to write ry(r) in terms of its time derivative p(r), which is a sum of delta 
functions at those times when phase slips occur. Because we take ry(r) = r]{T + (3), an equal number of "positive" 
(77 ^ 77 -I- 1) and "negative" (77 ^ 77 — 1) phase slips occur in each of the configurations summed over. In a configuration 
having n positive phase slips, we denote their times by r,, and the times of their negative counterparts by fj. Then 
we have: 



n 

PW = ^(<5(r-T,)-^(T-f,)). 



(54) 



i=l 



Without loss of generality, we can take ?7(r = 0) = and write r;(r) = JJ" dr' p{t'). In the Villain partition function, 
the sum over rj becomes a sum over p, which is easy to define precisely: 



E = E(^) J dn...dTn j dn...dfn. 



(55) 



p(r) n=0 
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Here, it was necessary to introduce an energy t„, which is the fugacity per unit imaginary time for phase-slip events. 
The factor of 2~^" is inserted for later convenience. Upon transforming to Fourier space and integrating out <p, we 
obtain a partition function only in terms of p: 



Zv = ^ exp 

Pir) 



1 + 



27rp(w„) 



(56) 



The duality transformation is completed by observing that Eq. (56) is equivalent to the following sine-Gordon 
theory: 



^dual 



VO exp ■ 



T 
2 



ln(^ 



ojl\0{uJn)? +U / dTCOs(27r6') 



(57) 



The equivalence is simple to establish: expanding in powers of ty and denoting the Gaussian part of the action by 
5q"*'[0], we obtain the partition function 



ual 



j.n f p 

Ti=0 CTi---cr„=±l 



P6lexp 



i=l 



(58) 



At this point it is helpful to note that 9{t) 0(r) + c is a symmetry of Sq^'^^, which means that only "neutral" 
configurations, i.e. those having X^^Uj = 0, contribute. This allows us to reorganize the above expression and write: 



^dual 



Wexp 



(59) 



Upon integrating out 6 we immediately recover Eq. (56). 

It is important to note that the term proportional to 1/J in the dual action is dominated by the logarithm at low 
frequency. Therefore, J will not not play a role in the low-temperature limit. However, it is important to keep this 
term in the formulation. This permits one to access the J <^ Ks and T ^ Tj regime discussed in Sec. I, where the 
resistance displays the truly activated behavior (i.e. having no logarithmic temperature-dependence) of Eq. (6). 

From these considerations, we see that the operator exp(27r«6'(T)) inserts a phase-slip event at time r. The ty cos{2tt9) 
term in Eq. (57) can therefore be interpreted as a vortex-hopping process. Furthermore, because phase and number 
are canonically conjugate variables, we can think of ^ as a charge difference between the films; in fact, it is precisely 
the charge imbalance Qimb between the two films {i.e. half the difference of the total charges of the films): 



(60) 



where Qi is the total charge on film i. Therefore, the operator for the voltage drop across the contact (from film 1 to 
film 2) is given by: 



V = -j^ (-f„cos(27r^)) = — — (f^cos(27r0)) =-^sin(27re). 



(61) 



B. Calculation of the resistance 



In order to find the resistance to lowest order in ty, we need to calculate the response of the voltage to a weak 
external current, using the dual partition function Eq. (57). This can be done by working in real time and replacing 
0^6 — It/2e in the vortex-hopping term [i.e. the tvcos{2n0) term]. Here, / is the externally-controlled current 
flowing across the contact from film 1 to film 2. The situation is now completely analogous to that of the standard 
derivation of the Kubo formula for a conductance, and we obtain for the resistance: 

(62) 

Here, G^{co) is the real-time retarded Green's function for the voltage: 

G'^io^) = r dte{t)e-^-'{[V{t),Vm)t.=o (63) 

J — CO 

= -i^ dte-'-*(e^W-e^(-*)), (64) 
2e Jo 



R{uj,T) = -Re -G)^(a.O 
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where 

F{t) = 4Tr^{e{t)e{0)-e^) (65) 



cos{xt) — l) coth — i sin(a;i) j 2- (^6) 

' " ■ UcK,/T^^J + \n{{uJc/xY-l)\ 



This correlator is obtained by the analytic continuation of its imaginary-time analog; see Appendix B3. After 
straightforward manipulations, the d.c. resistance can be written as 



R^T) = ^ / dtte^^K (67) 

We note that F{t) is analytic in the entire complex plane. 

Let us first consider the case of zero temperature, where we can write 

F{t) = WcK, / - 1) } 72- (68) 

We shall continue the resistance integral, Eq. (67), into the lower half-plane. If we let t = \t\e'^^ , where 7 G [— 7r,0], it 
is easy to show that, up to logarithmic corrections, 

Re[F(t)] ~-a(7)|t|, (69) 

for \t\ 00, where 0(7) > 0. Therefore, the integrand of Eq. (67) decays exponentially in this limit, and, by deforming 
the f-contour arbitrarily far into the lower half-plane, we see that R{T = 0) = 0. 

Next, we consider the case of nonzero temperature, for which it is convenient to write 

F{t) = ScK !^ i rei^(*+W2) , e-«(*+W2) _2cosh( — )1 - (70) 

Vo a;2sinh(f)L ^ ^ ^K^+[AcK,/n^J + \n{{oj,/x)^-l)]' 

Again we allow t to be complex and consider F{tr + iti) for ti G [— /3/2,0], and again up to logarithmic corrections 
we find that Re[F(t)] goes linearly to —00 for tr ±00. Then, instead of integrating along the real axis, we can 
integrate along the contour swept out hy t — i(i/2 for real t to obtain 

RiT) = l^j_Jtef^\ (71) 



where 



m = F{t-iP/2) = 16cKj -2 ■ UP^^ , u r. , 27^1 I ^2 .M^ ' ^^2) 

Jo X smh (i^-) + [AcKs/'k^J + In {{u)clxY - l)\ 

At low temperature, the resistance integral (71) can be evaluated by the saddle-point method. It is easy to see that 
/(O) — > —00 in the T — > limit and, furthermore, that /(O) > f{t) for all t (and T) - this means that f = is the 
global maximum of the integrand. We then have 



2j.2 /-oo -. 2^2 / 9 

R(T) « J^e/(°)/_^exp(lr(0)t^) = ^yz|:^exp(/(0)), (73) 

which is expected to be asymptotically exact as T — > 0. As shown in Appendix C 3, a conservative criterion for the 
validity of the saddle-point evaluation is that /(O) should be large and negative. This means that we should expect 
Eq. (73) to be a good approximation for the resistance whenever /(O) < —1. 

We then need to evaluate the low-temperature form of /(O) and /"(O) - this is done in Appendix C. The results 
are: 

/(0)«^£^ I =.Em (74) 
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and 

These relations are asymptotically exact as T ^ 0, where the dependence on J drops out. More precisely, the ratio of 
Eq. (74) to the exact value of /(O) goes to unity as T ^ 0. However, the difference of these two quantities does not 
go to zero; this is due to additive corrections to /(O) that are proportional to l/(Tln^(r)), which are sub-dominant 
but nonetheless diverge at zero temperature. These statements also hold for /"(O) and Eq. (75). 

We show in Appendix C 2 that Eqs. (74) and (75) are also asymptotically exact in the extreme good and poor 
contact limits. Precisely, in the good-contact limit we require that J ^ Kg and T <C Ks- In this case, wc recover the 
universal low-temperature forms for /(O) and /"(O), but for all T <^ Kg- In the poor-contact limit we require that 
J and Tj <^T ^J. We then have 

and 

/"(O) ~ -An^TJ. (77) 

Note that we must make the restriction T ^ Tj, because when T < Tj the resistance crosses over to the universal 
low-temperature result. 

The expressions Eq. (74) and Eq. (75) can thus be thought of as an interpolation between two limits where they are 
exact. Although this interpolation between good- and poor-contact limits is approximate, these formulas are valid to 
reasonable accuracy over a wide range of parameters. For example, using the parameters of the setup (with ground 
plane) discussed in Sec. VII, a straightforward numerical evaluation of /(O) shows that the formula Eq. (74) agrees 
to within 15% for T < 5K. However, if either the universal low-temperature form {J /Kg oo) or the poor-contact 
form Eq. (76) is used, there is only about 50% agreement at T w 5 K. 



(75) 



(76) 



VI. MECHANICAL ANALOGY 



The resistance formula Eq. (2) can be simply interpreted in terms of an analogy to an essentially classical system of 
springs at finite temperature. We focus here only on the exponential part of the resistance, which gives the dominant 
temperature dependence. We write this in the form 

" ^"K - il/cKg)Hl!j2T) + 2/nu ) 
In the limit of large Josephson coupling J (or at very low temperature), the resistance becomes 

Except for the logarithm, the resistance is activated in temperature. On the other hand, when J is small and the 
temperature is not too low, we have 

ii(T)ocexp(-^), 

which really is an activated form. 

The above crossover between two kinds of activated behavior suggests that we might be able to understand the 
resistance in terms of the thermally activated dynamics of an effective classical system. We shall obtain a description 
of such a system starting from the action for the phases $ j = ipi (0) at the contact - it will be instructive to retain the 
phases separately here, rather than focusing only on the phase difference (j). The action takes the form S = Sq+Sq+Sj, 
where 

-1 

\M^n)\^- (81) 



(78) 



(80) 



AcKsT 



E 
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15 



^2 ^1 



*2-*2 



FIG. 4: Depiction of the classical spring system defined by the action (84). The spring constants of the left and right springs 
are determined by the power-law phase fluctuations in the films [kf = 4cKs/it^ \n{LiJc/T)]. The central spring represents the 
Josephson coupling and has spring constant J. For a phase slip to occur, the system must cross over the energy barrier 
associated with stretching the coupled system a distance n. 



We wish to replace the slowly-varying logarithm by an a;„-independent constant. As LUn G 27rTZ, we have, very 
roughly. 



ln^^«21n-^. (82) 



The term Sj is the Josephson couphng, which we expand to quadratic order: 

r0 



(83) 



Furthermore, we imagine that the phase deep within each film is fixed to a classical value $i. 
With these approximations, we arrive at the action 



dr 







2cK - _ J 91 

^ -(($l-$l)2 + ($2_$2)2)+ . (84) 



7r2ln(w,/r)^' ' ^/ ■ V - ^' J ■ 2 



Note that we do not explicitly account for the 2tt periodicity; more precisely, if we had been more careful, the action 
would be independent of n if $2 = $i -I- 27rn. This action describes a simple mechanical picture of three springs 
coupled in series, as shown in Fig. 4. In order for a phase slip to occur across the contact, the system must tunnel 
over an energy barrier given by setting $2 ~ ^1 = This is equivalent to stretching the combined system of three 
springs by a distance tt, for which the excitation energy is found to be 

E = ^fceff(^)' = ^ (85) 

where we have used the rule for adding spring constants in series: 

In this picture, the resistance is proportional to the Boltzmann weight of the activation barrier exp(— i?/T), so we 
recover Eq. (78). We thus see that it is possible, at least roughly, to think of the low-temperature resistance as arising 
due to the thermal activation over an effective, temperature-dependent barrier. The properties of the barrier are 
determined by the long- wavelength physics of the superfluid films (and also their geometry, via the opening angle) . 
That one does not truly have simple activated dynamics is evident in the logarithmic temperature-dependence of the 
barrier height, which arises from the QLRO in the films. 



VII. RELATION TO EXPERIMENTS 



In this section we discuss the relevance of our theory to experiments on superconducting films. Although we consider 
only systems in the wedge geometry discussed explicitly in the paper, we expect that the basic features of our results 
will apply to a much broader class of Josephson tunneling experiments into and between 2d superconductors. Indeed, 
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we hope to stimulate more experimental and theoretical interest in such systems; some speculations on future directions 
such work might take arc contained in Sec. VIII. 

The most striking feature of the resistance formula [Eq. (2)] is the effective activation barrier Ea{T). In the 
universal limit, i.e. Kg/J<^ 1 or T ^ 0, this is set primarily by the superfluid stiffness in the films. The stiffness 
can be independently determined in thin film samples, via an inductance measurement,^^ so it should be possible 
to test this prediction rather thoroughly. Furthermore, Kg could potentially be varied in situ by tuning an in-plane 
magnetic field. The barrier also depends on ooc and J, which should really be considered fitting parameters. We 
can estimate uJc via iOc ~ w(|fc| = 27r/^), where ^ is the coherence length. In many cases, J can also be estimated; 
for example, consider a junction that is a constriction of length t and width w between two films each having bulk 
superfluid stiffness Kg. Then we estimate J « wKg/l. These estimates constrain the range of values of and J that 
may reasonably be used in fitting to experimental data. 

The most interesting regime is the universal limit [Eq. (5)], where the dependence on J, and hence on the details of 
the contact, disappears. In general, it may be difficult to access this regime by lowering the temperature; one needs 
\n{tujJc/'2.kBT) » AcKs/tt'^J, and the logarithm increases only slowly as T is decreased. Instead, it would probably 
be easier to fabricate contacts having large J /Kg; this could be achieved, for example, by fabricating short and wide 
constrictions between thin films. It would also be interesting to look for the purely activated behavior that should be 
present over a wide temperature-range when J / Ks is small. 

Another important consideration in fitting to Eq. (2) is that it may be better to use the more accurate form that can 
be obtained by numerically evaluating the appropriate integrals, especially when J Kg. As discussed in Sec. VB, 
Eq. (2) is obtained from 



where 



fit) = ['^"^'^ rfM cos{2ut/l3) - COSh(M) 1 

T Jo sinh(u) TT^+[AcKg/TT^J + \n{{0ioj2uY-l)]'' 

It is important to work with this more accurate form because the corrections to /(O) not contained in Eq. (2) are 
proportional to 1/Tln^(r); these arc dominated only rather weakly by the l/nn(r) term, and they also diverge at 
low temperature. We note, however, that the form (2) is already rather accurate over a wide range of temperatures. 
For example, using the parameters discussed below in the setup involving a ground plane, the approximate value of 
/(O) is within 15% of its exact value for all T < 5K. 

A crucial issue is to determine whether the system is indeed in a regime where our theory applies. To illustrate the 
important considerations, we consider a system that has already been fabricated and studied by Chu, Bollinger and 
Bczryadin,^^ in which two MoGc thin films are connected by a narrow constriction. The films have a characteristic 
linear dimension L and exhibit a bulk BKT transition at Tbkt ~ 5 K. For concreteness, we assume that Kg/kB ~ 10 K 
for T <C Tbkt, although in practice the stiffness should be determined by a direct measurement. Wc have opening 
angle 6 = n and Coulomb interactions {i.e. a = 2). The fabricated constrictions discussed in Ref. 22 have width 
w « 20 nm and length £ « 100 nm. The coherence length is about ^ w 7nm. With these parameters, we have 
huJc/kB « 2000 K and J/ks ~ 2K. 

We also consider a modification of the above setup, in which a superconducting ground plane is separated from 
the superconducting films by an insulating layer of thickness d = 100 nm. As discussed in Sec. IIIB, the interactions 
are screened down to a short-range form and, with the other parameters as above, wc have the superfluid velocity 
Vg = {l/K)-^l&iTe^Kgd « 3.8 x lO^cm/s. In this case, provided d > the cutoff is set by d rather than ^; that is, 
Wc = = 2-K/d). For our parameters, we have tuVc/kB ~ 2000K. 

In real systems, such as these, when does the resistance formula apply? First, the resistance should be in some 
sense small for our result to be a good approximation. More precisely, we require that Ea{T) / kBT > 1, so that terms 
of higher-order in ty are not important. For both parameter sets discussed above, this is satisfied for T < 5K. 

Second, the films must be in the thermodynamic limit. That is, they must be large enough that we can safely 
neglect finite-size effects, as well as the influence of the outside environment. This will be the case when plasmons 
undergo thcainal decoherence within the films much faster than they can travel across it. For films of characteristic 
linear dimension L, this is the case for T T^, where is a crossover temperature defined by 

kBTL = hui{k = 27r/L). (89) 



It should be noted that Tl decreases with increasing L. For acoustic plasmons, the definition of simply comes from 
a comparison of the thermal decoherence time It = 2Trh/kBT, and the time to travel across the system, tL = L/vg. 
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For plasmons having lo oc \fk dispersion, one has to replace by the group velocity of the fastest plasmons in the 
system, which is given by v™'"^^ ^ Ti^i^ ^ 2tt/L). This actually results in Eq. (89) but with an extra factor of 1/2 on 
the right-hand side. However, one should also impose the condition that the lowest-energy plasmons have significant 
thermal occupation, which holds above Tj, as defined in Eq. (89). In any case, it is more important simply to have 
an order of magnitude estimate for T/,. 

A separate finite-size effect is the Coulomb blockade phenomenon, where the transport is dominated by the charging 
energy Ec required to add a single Cooper pair to one of the films, and insulating behavior results. This physics is only 
important when T < Ec/ks- Noting that Ec oc for short-range interactions, and that Ec oc 1/L for Coulomb 
interactions, in both cases one has Tl ^ Ec/ks for large enough films. Therefore, it is enough to require that T ^ 
to eliminate finite-size effects, including the Coulomb blockade. 

With a superconducting ground plane as above, Tj, w 2 x 10~^ K for i = 1 cm. On the other hand, without the 
ground plane, w 1 K for L = 1 cm. There are two lessons to be learned from this. The first is that the ground 
plane is quite helpful in lowering T^. The second is that, even in with the ground plane in place, it is important 
to have rather large films for our theory to apply over a wide range of temperatures. Without the ground plane, it 
may be possible to achieve a modest reduction in Tl by sandwiching the system, above and below, by an insulating 
medium having a rather large dielectric constant e, as Tl oc e~^/^. 

VIII. DISCUSSION 

In this paper we have derived a general formula for the resistance of a point contact between two infinite super- 
conducting films. This quantity is a measure of the strength and nature of fluctuations in the superconductors. The 
low-temperature resistance is very nearly activated, with an activation barrier Ea determined only by the superfluid 
stiffness of the films and simple features of the geometry and interactions between Cooper pairs. Deviations from 
purely activated behavior enter via a weakly temperature-dependent reduction of the barrier height that depends only 
logarithmically on temperature. This result complements, and is intermediate between, the much-studied physics of 
3d bulk electrodes, where the Josephson effect occurs (R = 0), and Id superconducting wires, where Luttinger-liquid 
physics leads to a power-law tunneling resistance {R oc T"). In the 2d setting, the quasi- long-range order below Tbkt 
is not strong enough to enforce zero resistance at nonzero temperature. However, phase fluctuations are strongly 
suppressed, and the nearly activated resistance is quite close to a true Josephson effect. 

As discussed in detail in Sec. VII, it should be possible to test our results directly, via measurements of point- 
contact tunneling transport between superconducting thin films. More generally, we expect that the basic features of 
our results - most notably, the nearly-activated form of R{T) - will survive in a broader variety of systems involving 
tunneling into and between 2d superconductors and superfluids. It would be interesting to understand how this 
physics plays out in geometries more complicated than the one considered here, perhaps involving multiple contacts, 
multiple films, and other "circuit elements." 

In addition to enhancing our understanding of QLRO in usual superconducting films, superconducting point-contact 
measurements could perhaps be used as a probe of fluctuating off-diagonal order in poorly-understood systems, such as 
the high- Tc cuprates, and disordered fllms near a superconductor-insulator transition. Our work provides a baseline for 
such measurements, rooted in the well-understood physics of the superconducting state below Tbkt- Many interesting 
strongly-correlated and disordered materials possess gapless quasiparticle excitations - for example, disorder and 
field-induced quasi-particles in amorphous films or nodal quasiparticles in the cuprates — to which a point-contact 
measurement would be sensitive. In strongly disordered systems, the geometry could change in a nontrivial way, for 
instance if the QLRO in the film is sustained on a superconducting cluster near the percolation threshold. It would 
be very interesting to carry out analogous tunneling resistance calculations for the above cases. 

Even in the simple system considered here, several interesting issues remain to be addressed. One important piece 
of the physics not discussed here is the zero-temperature nonlinear I-V curve. This would give an additional measure 
of the off-diagonal order. It would also be useful to develop a better understanding how the phase-slip tunneling 
amplitude ty, which plays a crucial role in our analysis, depends on the more fundamental parameters of the system. 
Finally, we remark that a renormalization-group treatment has contributed substantially to the understand of the Id 
analog of this problem. So far, such an approach has not been successful here, but an analysis along these lines may 
be possible. 
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APPENDIX A: EFFECTIVE ACTION FOR THE POINT CONTACT 



In this appendix we derive the effective action for the phase difference across the point contact for the various 
geometries and interactions of interest. Most of these manipulations are completely straightforward; however, in the 
case of Goulomb interactions in the wedge geometry (Sec. A 4), we have found it necessary to modify the form of 
interaction to proceed analytically. As discussed in more detail in Sec. Ill, we believe that our modification does not 
affect the universal physics or the form of Eq. (13). 



1. Bilayer geometry, short-range interactions 

The imaginary-time action for the planes and contact is 5 = + ^2 + ^j, where 

/3 



* 



The fields ipi{r,T) are periodic in imaginary time, and have the following Fourier decomposition: 



/ (0^»(^'-")^" 



r — tLUriT 



(Al) 



(A2) 



<fi{k,u)n) = j dr j d'^r>fii{r,T).e-'''"'+' 



In Fourier space the action takes the form: 



d'^k /K 



)2 V 2 



(A3) 



Here, it is convenient to write the Josephson coupling using two auxiliary fields (^(r) and C('''): 



Jo 



dT{ iC 



•- {ipi{0,T) - ip2{0,T)) -JC0S{(I)) 



(A4) 



Integrating out ({t) identifies <P{t) with the phase difference across the point contact. 

It is now a simple matter to integrate out ipi and ( and obtain an action only in terms of 0: 



d^k 



(27r)2 fe2 + {oJn/Vs)^ 



Jo 



dr cos(/>(t). 



(A5) 



In Eq. (A5), we can easily carry out the k integral by imposing a short-distance cutoff for the modes in the planes, 
thus making the restriction |A;| < fcmax- Defining ojc = Vskmax, we obtain the point-contact effective action 



S = nKsTj2\k^n)\' 



In 



9 , 9 H — 1 



r0 

— J dr cos(/)(t). 
Jo 



(A6) 
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2. Bilayer geometry, Coulomb interactions 



In the case of Coulomb interactions, the action for one plane in terms of both the phase and density is: 



r 

Si = J drl J (fr{iif 



2\ , (2e)^ f cfrd^r'^"^^^'^'^'^'^ 



I 



r — r' 



(A7) 



Upon performing a Fourier transform in space but not imaginary time, this becomes: 



/f (Pk / 
dr J j^-^[iipi{k)ni{-k) 



KJk\ 



1 87re2 , 



2 Ifel 



(A8) 



Integrating out rn and taking the Fourier transform in the imaginary-time domain, we obtain an action for the phase 
variable 



S = 7E/7S.fS+A-,|.p)l*.P- 



T 

2^J (27r)2V87re2 
Proceeding as in the previous subsection, we find the effective action for (p: 



5 = ^^|^K)|^ 



d^k 1 



1 



(27r)2 |fc| (a;2/87re2)+i^,|fc| 



Jo 



dr cos0(r). 



(A9) 



(AlO) 



Performing the k integral (with the restriction |A;| < fcmax) and defining the high-frequency cutoff = Sne^Kskmax, 
the final result is: 



ln ""+"- 



-1 ./3 

— J / dr cos(/i(t). 
Jo 



(All) 



3. Wedge geometry, short-range interactions 

Here, we consider the wedge geometry in the case of short-range interactions, where the action for a single film has 
the same form as Eq. (Al). The only difference is that now, for each film, the position integration ranges only over a 
wedge; working with the plane polar coordinates (r, xp), this is simply the region Q < ij) < 9 and < r < oo. Rather 
than going to Fourier space, we shall transform to a basis of eigenfunctions x(^) of the Laplace operator satisfying 
the Neumann boundary conditions 



d^xir, "0 = 0) = d^x{r, ip ^0) =0. 



(A12) 



Letting — V^x = A^x, the eigenfunctions XXn{r) are defined in Eq. (21) and are labeled by a continuous variable 

A > as well as an integer n = 0, 1, 2, 

We can define the analog of the Fourier transform of the phase field ipf. 



^i{r,T) = / dA^e-*'^"^XAn(r)^.(A,n,u;„), 

w„ ■'^ n=0 

(Pi{X,n,uJn) = J dr J d^re"^''''x\n{r)'fii{r,T). 



(A13) 



In this basis the action & takes the form 



U,n Jo „=0 ^ ^ ^ 



Noting that 



f°° fx 
i{r = 0)= dX^-g>i{X,n = Q), 



(AM) 



(A15) 
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we can proceed as above to find the effective action for the phase difference across the contact: 



dX 



A 



X^ + iuJn/Vsy 



Jo 



dr cos 4>{t), 



(A16) 



where we have introduced the short-distance cutoff Amax- Defining Uc = VgAmax and performing the A integration, we 
obtain the final result: 



In 



2 1 2 1 ~ 1 



J I dr cos (j){T). 



(A17) 



Notice that this result is identical to that given in Eq. (A6), except for the factor of {9/2n), which accounts for the 

geometry of the wedge. 



4. Wedge geometry, Coulomb interactions 

Here, we replace the Coulomb interaction with the new potential V(r, r') discussed in Sec. III. The action for one 
plane, in terms of the phase and density, is 



/3 . 

Si = J drl J d^r{i^ini + ^iV^if)+^ J d^r d^r'V{r,r')niir)ni{r') 



Making the transformation of Eq. (A13), this action becomes 



S'i=T^^ / dX<uJri'f>i{\n,Un)ni{X,n,-Un) + 



47r6^ 

— \<f>i{X,n,uJn)\'^ + — ^|ni(A,n,a;n)P 



Upon integrating out the density, we obtain 

^^=2 EE / rfA(if.A2 + ^)|^,(A,n,a;„)|^ 
Proceeding as above, we find the effective action at the contact: 



^ = ^Ei^K)i' 



dX 



A 



A2 + {Xujl/8ne''K,) 



J dr cos(/>(t). 



(A18) 



(A19) 



(A20) 



(A21) 



The A integral is easily evaluated, and, defining = Sire^KsXmax, we thus obtain the final result for the effective 
action: 



5 = ^El^K)l' 



In 



2 I 2 "I ~ 1 



J dr C0S(/)(t). 
Jo 



(A22) 



APPENDIX B: ANALYZING THE IMAGINARY-TIME ACTION 



1. Analytic structure of the imaginary time action 

The action (13) for the phase difference across the contact has a rather peculiar form - the Gaussian part is 
logarithmic in frequency. In order to extract the correlation functions of the field (j){t) via analjrtic continuation, we 
need to understand the analytic properties of this logarithm. 

The function In can be broken into two parts: 



In 



2 I 2 
LU?, 



In h m 



(Bl) 
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FIG. 5: Analytic properties of In ^"'^ . The function has two branch cuts originating at lj„ = and terminating at lj„ — iiojc. 

" 2_| |2 

The values ±in resolve the value of the function near the branch cut, and are added to In ""^^j!^ to give the value of the 
function there. 



On the right hand side of Eq. (Bl), the first term has a branch cut stretching between a;„ = and ijJ„ = —itOc] 
similarly, the second term has a branch cut between a;„ = and aj„ — iujc- As one may expect, in order to be able 
to find correlation functions, we need to resolve the values of the logarithm functions in the vicinity of their branch 
cuts. The results are: 



ln-=^ 



In^ 



= In- 



LO—e-\-ix 



In 



uj— — €-\-ix 



UJn — — e — tX 



In 



= In- 



In- 



In- 



(B2a) 
(B2b) 
(B2c) 
(B2d) 



ujj, —e—tx 



where e > is an infinitesimal real number indicating which side of the branch cut we refer to, and < a; < lUc- The 

2 , 2 

properties of In a"*" in the complex plane are shown in Fig. 5. 



2. Phase correlation function 

In the weak-coupling limit, the equilibrium correlations of (f)(t) are related to the linear response conductance by a 
Kubo formula as described in Sec. IV. We are now in a position to find these correlations from the imaginary-time 
action, Eq. (13). 

For J = 0, the imaginary-time phase correlator is 



(</>(r)</)(0) - = r 



^ AcK. 



In 



!)■ 



(B3) 



As usual, we can rewrite the sum over Matsubara frequencies as an integral. This is done as follows: 



T 



y 



■In 



-1) = 



1 



AcK,, J 27r 1 - e-^""'^ 
C1+C2 



In- 



(B4) 
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FIG. 6: Contours of integration in the w^-plane for the Matsubara frequency sum Eq. (B4). Ci extends from e + ioo to e — ioo, 
and C2 ranges from — e — ioo to — e + ioo. Ci and C2 are closed by semicircles at infinity (not shown) in the right and left 
half-planes, respectively. The poles of (l — g-«"n/3^ ^ jjg qjj ^jjg j-gj^j axis, and are contained within Ci and C2, whereas the 
branch cuts of the logarithm function lie outside Ci and C2. 



In Eq. (B4) we introduce the contours Ci and C2 such that we include the poles of (l — e^"^"'^) but not the branch 
cuts of the logarithm function; these contours are shown in Fig. (6). The contours run infinitesimally to either side of 
the imaginary axis and are closed by semicircles at infinity. The only nonzero contributions arise near the imaginary 
axis for < uic, and we thus have 

(0(r)0(O)-02) = -^ 

where e = 0+. These integrals can be recast as integrals on the real axis by using the branch-cut resolution of 
Eq. (B2d): 

{Hr)m~^') - ^f-/— /— ^fln^^ + ^Tr) (B6) 

\ 

? idx e-""^ - \ ujI-x'^ \ 7 idx e-^^ - 1 /, tj^ _ ^.2 \ \ 

+ / a In -^—^ iTT - / 3- In -^^r h iTT . 

i 27r 1 - e^^f^ \ x^ y i 27r 1 - e-^/^ \^ ^2 ^ M 

/ 

Combining terms, we obtain 

(0(t)0(O) - <\?) ^ J dx (^(cosh(a;T) - 1) coth (^^^ - sinh(a;T)^ . (B7) 



The last step is to analytically continue this expression to real time, r — > it, which yields 

((^(i)</)(0) - (/.(O)^) = J dx (^{cos{xt) - 1) coth (^^^ ~ i sin(a:t)^ . (B8) 



We note that we have also obtained this result via the Hamiltonian formulation, working directly in real time. 



dlOr, 



27r 1 



-iuJnj3 



(B5) 



Dual correlation function 



In Sec. VA, we derived the dual form of the action at the point contact in terms of the charge-imbalance variable 
(r). In order to calculate the resistance, as discussed in Sec. VB, we need the correlation function {9{t)6{Q) — (P) t^=o. 
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This can be obtained by analj^ic continuation from its imaginary-time analog: 

J / e-*'^"^ - 1 



We can proceed exactly as we have just done for the phase correlation function. Rewriting the Matsubara sum as 
an integral over the contours Ci and C2 of Fig. 6, the result is: 

0K^ r" dx ( , , , , (i5x\ , \ 1 , ^ 

a Jo X y \2J y^2+(M^ + lni:I^) 

where the second equality is obtained by manipulations analogous to those of Appendix B 2. To analytically continue, 
we replace t ^ it and obtain the sought correlator: 

mm - 9') = ^ ^ ( ( cos(a;i) - 1) coth (^) - isin(xt) ) ^— ^ • (B12) 

APPENDIX C: EVALUATION OF RESISTANCE INTEGRALS 

1. Low-temperature form of /(O) and /"(O) 

Here, we derive the low-temperature forms of /(O) and /"(O), which go into the saddle-point evaluation of the 
tunneling resistance [see Eq. 73]. Recalling the definition of f{t) in Eq. (72), it is helpful to express it in the form 

_ d,cKs /■'^'^<=/2 ducos{2ut/l3) - cosh(M) 1 

~ T Jo V? sinh(u) 7r2+ [4c/<:s/7r2J + ln((/3a;c/2u)2- 1)]^' 



First consider /(O), which we write as 



/(0) = -^(/i(T) + /2(T)), (C2) 



where 

'■^ du cosh(M) — 1 



L(r) = / 
JO 



V? sinhu + [AcKJ-k'^J + In {{fiuc/2uY - l) 



(C3) 



and f2{T) is of the same form, but evaluated between the limits 1 and f3u}c/2. 

To evaluate fi{T), we first note that, to very good accuracy. In (^{(jwcl^uf' — l) « 2 ln(/3a;c/2u) for u G [0, 1]. We 
proceed by expanding [cosh(u) — 1]/ sinh(u) in powers of u. Keeping only the leading term (at small u), we have 

A(T) ' ^''^ 



8 Jo w (7r/2)2 + [2cKs/-k'^J + In {l3ujc/2u)] 
i|l_^arctan '^(\n{wc/2T) + 2cKs/t:'^ J^^Y (C4) 



When the argument of the arctangent is large. 



•^^(^^ ~ 8 ln(wc/2T) + 2cif,/7r2j ' (^^'' 
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This holds within 1%, as long as the argument of the arctangent is above about 4, which is always easily satisfied for 
reasonable parameters. The higher terms in the power-series expansion fall off as 1/ ln^(wc/r) or faster in the T ^ 
limit, so the T — > asymptotic behavior of fi{T) comes from Eq. (C5), and reads 

Now consider /2(T) - we will show that it is dominated by /i(T) as T — > 0. As [cosh(u) — 1]/ sinh(u) < 1 for all 

u G [l,/3a;c/2], we have 



r"^'/^ du l_ 

[AcKs/tt'^J + In {{(3ujj2uf - l)] ' 

We can break this integral into two more pieces, writing f2iT) < gA{T) + gsiT), where 

7r2 + [AcKs/tt^J + In {{(3u}j2uf - l)] 



h{T) < / ^ J. (C7) 



and gsiT) takes the same form but is integrated from to f}u}c/2. We take the low-temperature limit of 

gA+ gs, keeping e fixed. 

As long as e is not too large, we have 



1 J 

9a{T) < 2 / ^ (C9) 

TT^ +[4cK,/n^ J + In {{^,/ey-l)Y J 1 

1 



On the other hand, 

which is dominated by gA- Therefore, we have 



(CIO) 



gB{T)<— — ^J — , Cll 



In^K/T)' 

as claimed. Thus we may write 



cKs 1 
/(O) « -^^^^^^^^^^ ^ 2cK,/n-'J • 



We also have the low-temperature asymptotic form 

~rin(a;c/2r) 

Now we evaluate /"(O), by writing 



rK 

/(O) ~ (C14) 



/"(O) = -16c/f«T / . 

Jq smhw 7^2 + [4ci^,/7r2 J + In ((/3a;e/2w)2 - l)] ' 



-l&cK.Ti I du+ du 



sinhu + [4cK,/tt'^J + In {{f3u;j2u)'^ - l)] 
= -16ci^,r(/i(T) + /2(r)). (C15) 

The function /i(T) can be evaluated similarly to /i(T), by writing 1/sinhu = sinhu), and expanding 

u/ sinhu in powers of u. The leading term again gives the dominant T — > behavior, and proceeding as above we 
find 

\ln{uj,/2T)l2cK,/nU ' ^^^^^ 
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Furthermore, as 1/ sinh?i decays exponentially at large u, it is clear that f2{T) will go as 1/ hi^{u)c/T) at low temper- 
ature, and will be dominated by fi{T). Therefore, we have 

- -^--^-^ ln(../2T)!2c/../.V ' 

and, at low-temperature, 

2. Evaluation of /(O) and /"(O) for large and small J/K, 

The analysis above yielded the low-temperature asymptotic forms for /(O) and /"(O) valid for T ^ Tj, and hence 
the universal low-temperature limit of the resistance. Furthermore, we obtained approximate formulas for these 
quantities that reduce to the asymptotic forms for T <^Tj, and include some dependence on J away from this limit. 
Here, we show that these formulas are exact, not only at very low temperatures, but also in suitably defined good- 
and poor-contact limits. 

The good-contact limit is simply J ^ Kg. Taking J/ Kg — > oo, it is clear that all dependence on J simply drops 

out of /(i), and at low temperatures one can follow the analysis of the previous subsection to obtain the universal 
forms for /(O) and /"(O). The only requirement on the temperature in this case is that Kg/T ^ 1, so that Ea{T)/T 
is large, and terms of higher-order in ty do not contribute appreciably to the resistance. 

The poor-contact limit we must consider here is not merely the "naive" poor contact limit oi J <^ Kg. Because 
the resistance crosses over to the good-contact form at low enough temperature, we must also demand that T ^ Tj. 
With these conditions, the approximate forms of the previous subsection become: 

/(O) « (C19) 

and 

/"(O) w -Att^JT. (C20) 

These results are actually asymptotically exact for J ^ Kg and T ^ Tj. Furthermore, we should require that 
J/T ^ 1, so that the resistance is dominated by the contribution of leading-order in ty. 

To see this, let us revisit the analysis of /(O) given in Appendix C 1, and consider the fimctions fi(T) and fiiT). 
Again, fi{T) can be analyzed by expanding (cosh(M) — 1)/ sinhu in powers of u. Only the leading term in this 
expansion leads to a contribution to /i(T) of order J /Kg - this is exactly the integral we already evaluated to arrive 
at the approximate formula, Eq. (C13). All other terms in the power-series expansion give contributions of order 
(J/Kg)"^. Similarly, it is straightforward to show that f2{T) ~ (J/Kg)"^. Therefore, we can write: 

/(0)«-^ . (C21) 

2T i + ^ln(a;,/2T) 

Note that it is important to require that T l^- Tj; otherwise, {J /Kg) ln(u;c/2T) would not be small, and it would no 
longer be legitimate to neglect higher-order terms in J /Kg. However, for T ^ Tj we have: 

/(O) ^ (C22) 

Similarly, it can be shown that, in the same limit, 

/"(O) ~ -Att^JT. (C23) 



3. Validity of the saddle-point approximation 

We here derive a conservative criterion for the validity of the saddle-point evaluation of the resistance integral, 
Eq. (73). This criterion is simply that, in the power-series expansion of f{t), terms higher than 2nd order in t should 
be small when the 2nd order term is of order unity. Formally, we require that 

^|/(2")(0)|i^"«l, (C24) 
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where n> 1, and to is defined by 



2l/"(0)|io = 1- 



(C25) 



The higher-order terms alternate in sign, and wiU thus tend to cancel one another out, so we expect the saddle-point 
approximation to be even better than the results derived below suggest. 
First, we must estimate /^^"'(O), which we write in the form: 



/(2")(0) = 16cKs{-ir{2Tf''-H2n, 



where we have defined 



i2n 



L 



0i^c/2 ,,2(n-l) 



1 2 • 



(C26) 



(C27) 



sinhu + [4c/sr^/7r2 J + In ((/3a;c/2u)2 - l)] ' 
As the logarithm varies very slowly, and the integrand decays exponentially at large u, we can approximate as follows: 

1 r w^^"-^^ 



I2n 



7 

Jo 



7r2 + [4cK,/7r^J + 2 In (/3a;c/2u™ax)] sinhu 

1 (4"-2)(2n-2)!C(2n-l) 



7r2 + [4cK,/7r^J + 2 In {pLJc/2u^^^)y 



22n-l 



(C28) 



Here, Umax is the value of u at which m^^" sinhw is maximum; Wmax ~ 2(n — 1), to good accuracy. Furthermore, 
we can safely make the approximations 4" — 2 w 4" and C(2« — 1) ~ 1, so that 



(2n)! ^" ' n(2n - 1) [AcKs/tt'^J + 2 In (/3wc/4(n - 1))] ' 
Then, using Eqs. (C25) and (C17) to determine to, we have 



(C29) 



(2n)! 



|/('"nO)|io" 



8cK, 



'T[ln(wc/2r) + 2cKs/Tr^J] 



n{2n - 1)T tt^ + [icKs/n^J + 2 In {uJc/^T{n - 1))] ^ 
8(ln(cjc/2T) + 2cKs/tt^J) 



n(2n- 1) 



[AcK./n'^J + 2 In (a;c/4r(n - 1))] ' 



/(O) 



l-n 



(C30) 



The important part of this expression is the factor of |/(0)|^~" on the right hand side. We see that the condition 
Eq. (C24) for the validity of the saddle-point approximation is simply that |/(0)| is large or, equivalently, that the 
argument of the exponential in the resistance formula is large and negative. 



APPENDIX D: VARIATIONAL METHOD 



The main result of this paper is Eq. (2) for the contact resistance. This result, which was obtained by starting from 
the good-contact limit and implementing a perturbation theory in phase-slip events, is expected to apply for a large 
range of parameters that stretches all the way to J ~ T at low temperatures. 

As discussed in the Introduction and Sec. Ill, several pieces of evidence justify this expectation. Here, we provide 
further supporting evidence, obtained from a variational approach of the type first employed by Fcynman""^. 

In the context of Josephson or point-contact systems, the variational approach used here allows us to find a 
quadratic substitute to the cosine terms in the imaginary-time action. This method was used to find the phase diagram 
of the resistively shimted Josephson jimction,'^'^ and has since been used in several other contexts, e.g., quasiparticle 
dissipation''^ Although this method is an uncontrolled approximation, it may provide some information, if interpreted 
appropriately. 

The starting point of the variational calculation is the action in Eqs. (13) and (16), but written in terms of an 
Ansatz quadratic action plus a perturbation: 



S — So + Si, 



(Dl) 
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FIG. 7: Phase diagram resulting from the variational approach. A nonzero self-consistent solution to Eq. (D9) for the variational 
parameter D, exists for J > Jo{T), given by Eq. (DIO). The dashed line is the continuation of an high-T expression (only the 
explicit solution of the low-T regime is discussed in the text). 



with 



and 



In 



Si 



-I 



dr 



J(l-cos<^(r))--£)(^(r) 



(D2) 



(D3) 



The variational parameter is D, which is the coefficient of the ip'^ "mass" term in the quadratic Ansatz action Sq. 

To find the best 5*0 to describe the point-contact system, we follow the Feynman strategy by minimizing the free 
energy with respect to D: 



+ 



F = -lnZo + {Si)s^ = - ^[^] (-^o) 
A straightforward calculation yields 



jD[^]exp{-So) " 



(D4) 



(D5) 



where {■ ■ ■)so denotes the thermal average with respect to the action ^o- This leads to the following self- consistency 
condition for D: 



D = Jexp 

To evaluate the sum TY^g{uj)^ we follow exactly the same steps as in Appendix B, and thus obtain 



(D6) 



, , 1 f dx , f x3 



2 ■ 



(D7) 



The dominant contribution to the above integral comes from the region x < 2T, where the hyperbolic cotangent is 
large. We therefore make the approximations: cothx/3/2 « 2/{xP) and Incjc/x tt. This allows us to evaluate the 
integral approximately, 



2T 

4cK,T f dx 



1 



2cK,T 



1 



(D8) 
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and hence arrive at a transcendental self- consistency condition for D: 



Jexp 



TT^D 2£)ln^ + 



D. 



(D9) 



For all values of the parameters this condition has the trivial solution D = 0. It also has & D ^ solution for all 

temperatures, but only for J > Jo(T); i.e., the coupling must be stronger then a temperature-dependent critical value, 
Jo(T) (see Fig. 7). Because we are interested only in temperatures well below Tbkt, we focus on the low-temperature 
limit of the self-consistency relation, which is specified by Tln(a;c/2T) <C 2'kKs. In this limit, we find that JoiT) has 
the following form: 



^o(T) « -T, 



(DIO) 



where e is the base of the natural logarithm, not the electron charge. According to Eq. (D9), D jumps discontinuously 
at the critical line JoiT), from for J < Jo{T) to 



e 



(DU) 



for J > Jo(T). 

The appearance of a self-consistent mass term suggests that the Josephson coupling J should be treated as large, 
because the point-contact phase difference tends to localize in one of the wells of the Josephson potential. Therefore, 
the proper approach in the regime J > Jq{T) is to expand about the good-contact limit, where tunneling between 
minima of the Josephson potential is treated perturbatively. Indeed, the region where we find a self-consistent nonzero 
mass (Eq. DIO) essentially coincides with the regime of validity of the resistance formula (2), which is roughly given 

by 



T < Ea{T). 



(D12) 



This provides further justification that the expansion about the good-contact limit is valid at low temperature, even 
for small Josephson coupling. 
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